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Abstract 

Supersymmetry with R-parity violation (RPV) provides an inter- 
esting framework for naturally accommodating small neutrino masses. 
Within this framework, we discuss the lepton-flavour violating (LFV) 
processes ^ cy, fi — > eee, and /i ^ e conversion in nuclei. We 
make a detailed study of the observables related to LFV in different 
RPV models, and compare them to the expectations of R-conserving 
supersymmetry with heavy right-handed neutrinos. We show that the 
predictions are vastly different and uniquely characterise each model, 
thus providing a powerful framework for experimentally distinguishing 
between different theories of LFV. Besides the obvious possibility of 
amplified tree-level generation of /U ^ eee and /x ^ e conversion in 
nuclei, we find that even in the case where these processes arise at the 
one-loop level, their rates are comparable to that of // — > 67, in clear 
contrast to the predictions of R-conserving models. We conclude that 
in order to distinguish between the different models, such a combined 
study of all the LFV processes is necessary, and that measuring P-odd 
asymmetries in polarised fi — > eee can play a decisive role. We also 
comment on the intriguing possibility of RPV models yielding a large 
T-odd asymmetry in the decay of polarised fj, eee. 



1 Introduction 



Recently, neutrino oscillation experiments 0, 0, 0] have provided very strong 
evidence for non-zero, yet tiny, neutrino masses. In order to accommodate 
such small masses, it is widely believed that new physics beyond the Standard 
Model (SM) is required. One of the simplest and most elegant mechanisms 
for generating a small neutrino mass is to introduce extra standard model 
singlets to the SM Lagrangian, and allow them to acquire a very large 
Majorana mass (this is the well known seesaw mechanism [^). There are 
many important phenomenological consequences of neutrino masses. One 
of them is that individual lepton-flavour numbers are not conserved, which 
implies that SM forbidden processes such as fi —>■ ej may occur. However, 
given the size of the neutrino masses, the rates for charged lepton flavour 
violating (LFV) phenomena are extremely small in the SM plus massive 
neutrinos 

There are other hints for physics beyond the SM, including the gauge hi- 
erarchy problem. Low-energy supersymmetry (SUSY) is one of the preferred 
candidates for beyond the SM physics which solves the hierarchy problem. 
SUSY models can easily accomodate the seesaw mechanism, and SUSY even 
helps in the sense that it stabilises the (very heavy) Majorana mass of the 
right-handed neutrino. Furthermore, in such a framework, LFV processes in 
the charged lepton sector such as /i ^ 67, /i eee, and fi ^ e conversion in 
nuclei are potentially amplified, as has been previously discussed |^, ||, |[ , 
and the rates for such processes can be within the reach of future experiments. 
The reason for this is that while in the SM plus massive neutrinos the 
amplitudes for LFV violation are proportional to the neutrino masses {i.e., 
suppressed by the very large right-handed neutrino masses, in the case of 
the seesaw mechanism), in SUSY models these processes are only suppressed 
by inverse powers of the supersymmetry breaking scale, which is at most 
0(1) TeV. 

Another SM extension which naturally accommodates non-zero neutrino 
masses is SUSY with R-parity violation (RPV). R-parity is usually imposed 
as a global symmetry of the minimal supersymmetric version of the SM 
(MSSM) in order to prevent an unacceptably large rate for proton decay. 
However, this proves to be somewhat of an overkill, since R-parity conser- 
vation implies both baryon number and lepton number conservation, while 
to stop proton decay only one or the other needs to be exactly conserved. 
In light of the evidence for neutrino masses, which can potentially be Ma- 
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jorana particles and therefore violate lepton number, one may, instead, take 
advantage of RPV operators to generate small neutrino masses. 

In this paper, we consider SUSY models with RPV but with baryon 
parity (in order to satisfy the current experimental upper limits on the proton 
lifetime These models naturally generate small Majorana neutrino 

masses, if the RPV couplings are small [|13|, In such RPV models, 

"large" LEV in the charged lepton sector is also generically expected. Indeed, 
as has been pointed out in the literature ^7\, ^ |22|, the most 

stringent limits on certain products of RPV couplings come from the present 
experimental bounds on charged LFV processes. Therefore it is important 
to understand some of the general features of LFV in models with RPV. 

It is interesting to consider how searches for LFV at low energy ex- 
periments compare to those at colliders. For instance, the similtaneous 
presence of R-violating operators that couple both to e — g and to fi — q 
(r — q) pairs, would lead to yU+jet (r+jet) final states at HERA p3[. It 



turns out that for e ^ r transitions, the high energy experimental probes 
provide the strongest bounds, while for e <-> /i transitions stopped muon 
experiments provide, by far, the most stringent bounds. Finally, the strongest 
bound for r ^ /i transitions comes from r — >■ /i7 searches at CLEO |2| 



[Brij fi'-f) < 1.1 X 10~^) which is less restrictive. In the near future, the 
experimental sensitivity to some rare muons processes is going to improve by 
two to three orders of magnitude, while a similar improvement is not expected 
for other LFV processes. For this reason, we will focus on processes with 
stopped muons, which not only provide the stringest quantitative bounds on 
LFV today, but which will be significantly probed in the near future. 

In this paper, we discuss the LFV processes fi^ /i"*" — > e~^e~e'^, 

and yU~ —>■ e" conversion in the case of models with trilinear RPV. In Sec. 2, 
we briefly introduce the SUSY models with trilinear RPV which will be 
considered here. In Sec. 3, we present the formalism for computing branching 
ratios and asymmetries of the relevant LFV processes. In Sec. 4, we consider 
LFV processes in some representative cases, including those in which the 
branching ratio for /i+ — e'*'7 is much smaller than the branching ratio for 

^In cosmology, large RPV Yukawa couplings may erase a pre-existing baryon 
asymmetry . Here we do not consider such constraints since they are model-dependent 
and can be evaded in several baryogenesis scenarios . 

mechanism which explains why RPV couplings are small is required. This can be 
achieved, for example, by imposing flavour symmetries which relate the lepton and baryon 
number violating Yukawa couplings to those that generate fermion masses . 
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/i"*" — i> e~^e~e~^ and/or the rate for ^~ e~ conversion in nuclei, which can 
be generated at the tree-level. Even if all LFV processes occur at the one- 
loop level, the rates for all three processes considered here are comparable. 
These features are completely different from the predictions of other neutrino 
mass generating SUSY frameworks, such as the MSSM with right-handed 
neutrinos 0. In the latter, the branching ratio for fi^ e+7 is much 
larger than that for —>■ e^e~e^ and the rate for yU~ conversion in 

nuclei, even though all processes are also generated at the one-loop level. 
We also show that P-odd asymmetries in the /i"*" — >■ e^e~e^ process (which 
require polarised muons in order to be measured) are very useful in order to 
distinguish different models. Sec. 5 contains our conclusions. In Appendix A 
we provide explicit expressions for the LFV vertices in the case of models with 
RPV, while in Appendix B we discuss the current bounds on certain pairs of 
RPV couplings from LFV processes and comment on neutrino masses. 



2 SUSY models with trilinear R-parity vio- 
lation 

Here, we briefly introduce the SUSY models with RPV which will be dis- 
cussed in the upcoming sections. If R-parity conservation is not postulated, 
in addition to ordinary Yukawa interactions, the following terms are allowed 
in the MSSM superpotential: 

^RPV = -^LiLjEf, + X'ijj^LiQjDk + X'ljJJiDjDk + fi'^LiHu, (2.1) 

where Lj, Ei, Qi, Ui, Di, and Hu denote the left-handed doublet lepton, 
right-handed lepton, left-handed doublet quark, right-handed up-type quark, 
right-handed down-type quark, and "up-type" Higgs superfields, respectively. 
The indices i,j and k range from 1 to 3 for different quark/lepton flavours. 
Throughout this paper, in order to forbid rapid proton decay, we impose 



baryon parity |10], so all A" couplings are zero. We also make the simplifying 
assumption that all fi' also vanish.^ In light of these assumptions, the 
superpotential above yields the following Lagrangian: 

■'Even if ^[ were non-zero, their contributions to LFV processes would be, in general, 
negligible because of neutrino mass constraints llj, 02] . 
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+ Kjk^KM {yLidLad*jik + dRk^Lidba + dji^dia^Li) 

- Kjk {i^'jeLid*Rk + dRk^LiULj + dRkULjeu) + h.c, (2.2) 

where f {f = v-, e, d, and u) denotes fermions and / sfermions, and the 
index (i?, L) indicates the field's chirahty. We assume that the RPV Yukawa 
couphngs above {Xijk and A'j^^) are the only source of LFV. In what follows, 
Eq. ( p.2|) is what is referred to by "RPV model." 

3 Branching ratios and asymmetries for the 
LFV processes 

In this section, we present complete expressions for the branching ratios for 
the LFV processes /i^ — > e"^7, e~^e~e~^, and fi~ —>■ e~ conversion in 

nuclei, for the P-odd asymmetry in /i+ — > e"'"7, and for the P-odd and T-odd 
asymmetries in — >■ e^e~e^ . 



3.1 /i+ e+7 

The process /i"*" e"*"7*^*) is generated by photon penguin diagrams (see the 
penguin diagrams in Figs. |1|, 0, and ^). The amplitude for this process can 
be written as follows: 



+m^iao,pq^{A^PL + A^Pr)] v,{p - g)(3.1) 

where u^(e) and e are the antimuon (positron) and photon wave functions, 
and p and q are the antimuon and photon momenta, respectively. P^ and 
Pr are chirality projection operators: Pl = (1 — 75)/2, and Pr = (1 + 75)/2, 
while (Tq,/3 = (z/2)[7q,, •y^]. The effective couplings Ai'^ come from off-shell 
photon diagrams (g^ ^ 0), which only contribute to fi^ e^e~e^ and 
— *• e conversion in nuclei. On the other hand, the couplings Ag'^ arise 
from the on-shell photon diagrams (g^ = 0), which induce /i+ — >■ e"^7 as well 
as /i"^ — i> e^e~e^ and /i~ ^ e~ conversion in nuclei. Explicit expressions for 
Ai'2 in models with RPV are presented in Appendix A. 



In the yU+ — > e"'"7 decay, it has been argued |^ that a nonzero muon 
polarisation is useful not only to suppress background processes, but also to 
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distinguish between yU"*" e^7 and /i"*" 6^7. The differential branching 
ratio for yU"*" — > e'^'y is given by 

= Mi:^£:2){i+^,P„,9}. (3.2) 

a cos y 2 

where P is the muon polarisation and 6 is the angle between the positron mo- 
mentum and the polarisation direction. Here, the branching ratio Br(yU"'" — *■ 
e"'"7) and the P-odd asymmetry Ap are 

Br(/.+ -e+7) = ^ (l^^P + I^^P) , (3.3) 

where is the Fermi constant, and a is the fine-structure constant. 
3.2 Polarised /i+ ^ e+e+e^ 

In the RPV models, some of the LLE couplings {Xijk) generate /i"^ e"'"e~e+ 
at tree- level (Fig. p, while the photon penguin vertices Ai'2 also contribute.^ 
The amplitude for /i"*" — > e^e~e"'" is 

+ B^Vf, (p) PR-ff,Ve {p2)Ue {p3 )PLl^Ve{pi) 

+4™^(p) ^{A[Pl + P^)7^ + m^i^iA^Pn + A^Pi)| ^^e(P2) 

XMe(P3)7''^^e(Pl) " {Pl ^ P2) , (3.5) 

where the explicit expressions for the tree-level vertices B^'^^^ in models with 
RPV are given in Appendix A. 

When the muon is polarised, two P-odd and one T-odd asymmetry can 
be defined |^ Using the notation introduced by Okada et al. [^, the 



2;-axis is taken to be the direction of the electron momentum and the {z x x)- 
plane is taken to be the decay plane. The positron with the largest energy 



^There is also a Z-penguin contribution. However, its contribution is suppressed by 
mj/m^ where m/ is the typical fermion mass in the process. Therefore we simply neglect 
it. In order to be consistent, we will not study processes where top-quarks are involved. 
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is denoted as positron 1 and the other as positron 2. The x-coordinate is 
defined as {pi)x > where pi is the momentum of positron 1. It is in this 
coordinate system that the direction of the muon polarisation P, used below, 



is defined. (For details, see ||2^.) Finally, the P-odd and T-odd asymmetries 
are defined as follows: 



A 



Pi 



A 



P2 



An 



N{P, > 0) - N{P, < 0) 
iV(P, >0) + iV(P, <0)' 

{0.61(Ci - C2) - 0.12(^3 - C4) + 5.6(^5 - C,) 

-4.7(^7 - Cs) + 2.5(C9 - Cio)} , (3.6) 
iV(P, > 0) - iV(P, < 0) 
iV(P, > 0) + iV(P, < 0) ' 

^{0.1(^3-^4) + 10(^5-^6) 

+2.0(C7 - Cs) - 1.6(^9 - Cio)} , (3.7) 

N{Py > 0) - N{Py < 0) 
N{Py > 0) + N{Py < 0) ' 

= ^l^{2.0Cn-1.6Ci2}, (3.8) 

where the muons are assumed to be 100% polarised, and N{Pi > (<)0) 
denotes the number of events with a positive (negative) Pj component for 
the muon polarisation. Here an energy cutoff for positron 1 is introduced 
(£^1 < (m^/2)(l — 6)) and henceforth we will consider 6 = 0.02, following 
Okada et al. |2^. This choice is made in order to optimise the T-odd 
asymmetry. Of course, one can obtain more information concerning the Cj 
coefficients, including the CP-odd terms Cu and Cu, (see definition of Ci in 
what follows) by analysing the Dalitz plot of the yU+ — > e~^e~e'^ decay. Ci 
(2 = 1 — 12) are functions of the effective couplings Ai'^ and B^'^: 

C. = = '^\Ak\\ (3.9) 

C3 = ^\B'' + A-naA^W = -^\B'^ + A-naA^W (3.10) 

oLx p o(J p 



c, = ^1^2 r, ^6 = ^|A^i^ (3.11) 
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^10 = -i^Re{yl^(i?«*+47raylf)}, (3.14) 



C, = -^Re{A^A[*), Q = -^RelA^Mf*), (3.12) 

Cg = -^Re{A^(S^* + 47ra^f*)}, (3.13) 

4G^ 

Cn = ^Im{87ra(^fAf + A^^f)}, (3.15) 

Ci2 = -^Im (5^* + Ai^aAt) + ^^(5^* + 47ra<*)} . (3.16) 
4G^ ^ ' 

The branching ratio for 5 = 0.02 is 

Br(5 = 0.02) = 1.8(Ci + C2)+0.96(C3 + C4)+88(C5 + C6) 

+U{Cr + Cs) + S{Cg + Cio). (3.17) 

The branching ratio for /i'^ — > e'^e'e'^ for 5 = is given by 

log^- — Wa + Ce) 



+ 16(C7 + C8)+8(C9 + Cio) (3.18) 



8^1 



+1 (log ^ - x) (1^2^^ + l^^n - 4Re(^[A«* + AM*)| 
+87rQ;Re {A^B^* + A^B^* - 2{A^B^* + A^B^*)}] (3.19) 

3.3 //~ — ^ e~ conversion in nuclei 

Similarly to — > e'^e~e"^, not only photon penguin diagrams but also tree- 
level diagram induced by some of the LQD Yukawa couplings (A^^-;;.) can 
generate ii~ — > e" conversion in nuclei. The amplitude is given by 

- Anaue I^iMi*Pl + A^*Pn) + m,i^{A^*Pn + A^*Pr^)^u, 

X Qq'^qlu'^q^ (3-20) 

q=u,d 



where the complete expressions for the tree-level contributions D^''^ in the 
case of RPV models are presented in Appendix A. The yU~ e" conversion 
rate is 



.^^^^''""^ [U^a^^ - A^\' (3.21 



+ 



IGtt^Z r(/i capture) 



Ai 



where r(yU capture) is the muon capture rate in the nucleus of interest p8[, Z 
and N are the proton and neutron numbers, respectively, F{q) is the nuclear 
form factor as a function of the momentum transfer and Z(.ff is the nuclear 
effective charge In some of the most commonly used nuclei, IgTi and 



27 
13 



Al, these nuclear parameters are given by 



r(yU capture) = 
Z = 22, Z,fj = 17.61, \F{q 



2.590 X 10^ s~ 

2 . 



1.7 X 10-1^ GeV 



and 



r(/i capture) : 

Z = 13, Zeff 



0.7054 X 10^ s"^ 



0.535 foriTi, 



4.6 X 10"^^ GeV 



(3.22) 



11.62, \F{q'' = -ml)\ 



0.64 for ^^Al. 



(3.23) 



4 LFV in representative cases 

The most severe constraints on some particular products of trilinear RPV 
couplings come from the present experimental upper limits on the branch- 
ing ratios of the LFV processes discussed in the previous section (see Ap- 
pendix B) |T3, ig, |T|, ^, ^ Therefore, searches for LFV in muon 
processes are particularly sensitive to models with RPV. Generically, it is 
very hard to make definite predictions for the branching ratios of the LFV 
processes since the number of new Yukawa couplings {Xijk, Kjk) is too large. 
Here we consider, instead, different cases where only a small number of RPV 
couplings is significant for LFV. This is done not only to simplify the problem 
at hand, but also to identify features of LFV which are not only different from 
those in the "traditional" models, such as the MSSM with heavy right-handed 
neutrinos discussed in P, |3, H, but which can also be used to characterise 
the different cases themselves. This dominance of specific RPV couplings is 



also a consequence of certain flavour models [16 
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A231 ^ ^Ai3i A231 ^ ^ ^^Ai3i L 

• • M — J — ^ — — " — 1 _ — r-^ 




Figure 1: Lowest order Feynman diagrams for lepton Havom violating 
processes induced by A131A231 couplings (see Eq. 



4.1 



e^e e 



induced at tree-level 



First, we consider a model in which only the Yukawa couplings A 131 and A231 
are non-zero. In this case, /i"*" — > e~^e~e~^ is generated at the tree-level, while 
the other LFV processes (/x^ e^7 and /i~ e~ conversion in nuclei) 
are induced via photon penguin diagrams at the one-loop level, as shown in 
Fig-i 

The effective vertices are given by 



A131A231 






2ml ' 






A131A231 


(- 




967r^m? 


2ml^ 


A131A23I / 


8 

~3 ~ 


2 log- 



mf^ 



+ 



mt 



25{ml/q^ 



(4.1) 
(4.2) 
(4.3) 



Here we assume, without loss of generality, that the RPV couplings are real. 
The function 5 is presented in Appendix A. In /i~ — > e~ conversion, we 
assume the momentum of the virtual photon to be = — in order to 
compute 5{m1/q^) in , while in the case of ji^ — > e^e~e^ , we simply set 
= (5 = 0), since the tree-level contribution is much larger than 
the contribution from 5{m^^/q^). The ratios of branching ratios, Br(yU^ — >■ 
e"'"7)/Br(yU+ — > e+e^e+) and R(/i^ e~ in nuclei)/Br(/i+ — > e+e^e+) do 
not depend on the R-parity violating couplings A131A231, and only depend on 
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the SUSY mass spectrum through m?^ and 



Brdu"*" ^ e"'"7) 
Br(/i+ — i> e+e~e+) 
R(/i- ^ e" in Ti (Al)) 
Br(/i+ — s> e+e~e+) 



4 X 10 




(4.5) 



the second of the equal signs being vahd for m,>^ = mg^ = 100 GeV. Here 
(3 = H-(one-loop contr.)/ (tree- level contr.) in the /i"*" — > e"^e~e+ process, 
which is close to unity (for example, P = 0.98 for = mg^ = 100 GeV). 
Since the fi~^ —>■ e+e'e^ process is generated at tree level, its branching ratio 
is much larger than that of the other LFV processes, as expected. If such a 
scenario were realized in nature, the /i+ e+e^e+ process would dominate 
over all the other channels, i.e., it is very likely that if nature realizes this 
particular scenario, n'^ — >■ e~^e~e~^ is within experimental reach while fi~^ — >■ 
e"'"7 is orders of magnitude below any foreseeable future experiment. 

Another interesting feature of Eqs. (f4.4|J4.5| ) is that, because of an ultravi- 
olet log-enhancement of the off-shell photon penguin contribution (Af ) ||2^, 
pO[ , the fi~ e~ conversion rates are significantly larger than the branching 
ratio of fi^ -—>■ e^'j. 

It is important to emphasise that the ratios of branching ratios of the 
different processes are very different from those in the different neutrino-mass 
models. For example, in the MSSM with heavy right-handed neutrinos (and 
R-parity conservation) 0, the following relations are approximately satisfied, 
because the on-shell photon penguin contribution A2 tends to dominate over 
all others: 



■^l 11 



Br(/i+ ^ e+e-e+) ^ {log — - ^ 

R(/i- ^ e- in Ti) a^Z^ffZ\F{q)\''mlGl 



1.6 X 10^ (4.6) 
0.92 (4.7) 



Br(/i+ ^ e+e-e+) 4^2 (i^g ^ - ^) r{fi capture) 

Another interesting feature of the case at hand is that in the /i"*" — > e"'"e~e+ 
process we obtain the following P-odd asymmetries (Eqs. ( |3.6| - |3.7| )): 
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-1.3, (4.10) 



A 



Pi 



Ap, 

since the tree-level contribution (C4) is dominant. The key feature here 
is that the two different P-odd asymmetries have opposite sign; Ap-^/Ap^ ~ 
— 1.3. More generally, this feature is present whenever the effective vertices 
j^L,R g^j^g dominant. In Table we list results of other similar examples. 

The situation is clearly different from the MSSM with heavy right-handed 
neutrinos, where the on-shell photon contributions A2 (C5) are dominant 
(this case is also listed in Table |1|, in order to facilitate comparisons): 

. 3(5.6C5) . 3(10^5) 



Ap, 



0.6. (4.11) 



Therefore, a measurement of the (sign of the) ratio of P-odd asymmetries 
in /i"*" —y e~^e~e'^ can clearly separate these two models (-Bl,_r ^ v4f''^ versus 

Another useful observable which may be measured in the case one has 
access to polarised muon decays is (Eq. (|3.4] )). In RPV models, Ap can 
have different values (see Table |l|), while in other SUSY extensions of the 
SM, either /i+ — > e^7 or /i+ e^7 is forbidden. Some examples include R- 
parity conserving SUSY with right-handed handed neutrinos (see Table P, 
SU (5) and 5*0(10) SUSY grand unified theories, and other MSSM extensions 



4.2 All processes induced at one-loop level 

Here we consider a different representative case, in which all of /x"*" — > e"'"7, 
/i"*" e~^e~e'^, and /x~ — > e~ conversion in nuclei are induced at the one-loop 
level (at the lowest order in RPV couplings) through the photon penguin 
diagram (Fig. H). Suppose, as an example, that only the couplings A132 and 
A232 are nonzero (again we assume that both of them are real, without loss 
of generality). The effective vertices for the LFV processes are 
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Table 1: The ratios of branching ratios Bi{ fi^ — >■ e~''7)/Br(/i+ — > e+e^e"*") and 
R(/i~ — s> e~ in Ti)/Br(/i"'" — >• e"'"e~e+), P-odd asymmetries Ap for /i+ — > e"'"7, 
and y4p2 for /i"*" — e"'"e~e"'" are shown when the hsted pair of Yukawa 
couphngs is dominant. Case (1), (2), (3) refers to the representative classes 
of models discussed in Sees. 4.1, 4.2, and 4,3, respectively. Here, we assume 
m- = 100 GeV and no mixing in the charged slepton mass matrix, and 
rriq = 300 GeV. We also show a typical result obtained for the MSSM with 
heavy right-handed neutrinos and R-parity conservation 





Br(/i-+3e) 


R{H^e in Ti) 
Br(/j^3e) 


Ap 




Ap, 


ApJAp, 


Case (1) 














A131A231 


1 X 10-^ 


2 X 10-3 


-100% 


+19% 


-15% 


-1.3 


^121-^122 


8 X 10-^ 


7 X 10-3 


+100% 


-19% 


+15% 


-1.3 


-^131-^132 


8 X 10-^ 


5 X 10-3 


+100% 


-19% 


+15% 


-1.3 


Case (2) 














A132A232 


1.2 


18 


-100% 


-25% 


-5% 


5.6 


^133-^233 


3.7 


18 


-100% 


-25% 


-4% 


6.2 


^231-^232 


3.6 


18 


+ 100% 


+25% 


+4% 


6.2 


^122^222 


1.4 


18 


-100% 


-25% 


-4% 


5.7 


^123^223 


2.2 


18 


-100% 


-25% 


-4% 


5.9 


Case (3) 














^lll'^211 


0.4 


3 X 10^ 


-100% 


-26% 


-5% 


5.4 


^112^^212 


0.5 


8 X 10^ 


-100% 


-26% 


-5% 


5.4 


^113^213 


0.7 


1 X 10^ 


-100% 


-26% 


-5% 


5.5 


^^121^221 


1.1 


2 X 10^ 


-100% 


-26% 


-5% 


5.6 


MSSM with Up 


1.6 X 10^ 


0.92 


-100% 


10% 


17% 


0.6 
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A232 ^ ^Ai32 A232 ^ ^Ai32 

• • M _ • _ J — ^ — — " 1 _ • _ I — " 




e (g) e (g) 



Figure 2: Lowest order Feynman diagrams for lepton flavour violating 
processes induced by A132A232 couplings (see Eq. 



(4.12) 



= 


A132A232 


(■ 






2m? 


A[ = 


A132A232 / 


8 
~3 


-21og- 



S^^-W)). (4.13) 

Again, we set = — for /i~ ^ e" conversion, and = for yU"*" — >■ 
e"'"e~e"'".0 The ratios of branching ratios Br(/i"'" —>■ e+7)/Br(/i+ ^ e'^e~e"^) 
and R(/i~ — e~ in nuclei)/Br(/i+ — >■ e"'"e~e^) are independent on the choice 
of A132A232: 

{-] -i. 

Br('u+ ^ e+7) , I 2^ 

- 3.2 X 10^ ^ 



Br(/x+ ^ e+e^e+) 



1.2, (4.14) 



R(/.- - in Ti (Al)) ^ 19.5 (n.5) ^ ^ ^ ^ ^ ^ 



1? m2 



2 



Br(/i+ — i> e+e~e+) 



3 21og^ + £t)%' 



8 



= 18(11), (4.15) 

'^Since the log-term is much larger than the 6 term for < (7^ < to^ in Eq. ( [4.13[ ), the 
result does not depend significantly on the choice of q in the e^e~ process. 
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where the second of the equal signs holds for rrip^ = rrijXj^ = 100 GeV. Here 
7 is a function of the SUSY mass spectrum, but it is of order unity. 

l\A^\'(tog§-ii)-4TtMAiA^) 
■'-'^ • \Ai\' ■ 

As an example, 7 = 1.09 for m/>^ = m^^ = 100 GeV. 

Because of the ultraviolet log-enhancement of the off-shell photon penguin 
diagram (Af ) in Eq. (|4.13|) , the event rates for the /i+ — > e~^e~e~^ and fi~ — >• e~ 
conversion in nuclei can be as large as the branching ratio for the /z^ — >■ e^7 
process, even though they are higher order processes in QED.0 Fig. |^ depicts 
the dependence on the slepton masses of these ratios of branching ratios. In 
the case of /i+ e+7, a cancellation between the two different diagrams 
(sneutrino and smuon loops) can occur, such that its branching ratio can 
be much smaller than that of the other processes. On the other hand, the 
numerical value of the ratio R(/i~ — > e~ in nuclei)/Br(yU"*" — > e+e~e+) is 
stable in a large region of the parameter space. All the LFV processes are 
equally relevant in this model. Again, we stress that these ratios of the 
branching ratios are very different in more "traditional" cases, such as in the 
MSSM with heavy right-handed neutrinos (see Eqs. ([4.6| , |4.7|) ). 



Since the off-shell photon diagram is dominant in the — * e+e^e+ 
process, C2 (= C4) in Eqs.( p.9| , p.l0| ) is much larger than the other Ci {i 7^ 
2,4). In this case, the P-odd asymmetries behave as follows: 

^ 3(-0.61C, + 0.12c,) ^ _ 

- 2(1.8^2 + 0.96^4) ^° ^ ' 



A 



0.19. (4.19) 



These relations (Eqs. (|0!^ , [4. 15| , ^. 191 ) ) are a typical feature of models in which 



the off-shell photon diagram is the dominant contribution to the — *■ 
e~e^e~ process. The results of other similar examples are also listed in 
Table |l|. 

^In the case of /i+ e+e^e+, t here is also an infrared log-enhancement to the 



branching ratio, as can be seen in Eq. ( 3.19 ) 
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Figure 3: Contours of constant Br{fi'^ — > e"'"7) /Br(/i"'" — > e"'"e~e+) ftopj, and 
R{fj,~ — >■ e~ in Ti) /Br{fi~^ — >• e+e~e"'") (bottom) in the (rrijXj^ x mj>^j plane, 
assuming that only the product of LLE couplings A132A232 is non-zero (see 
Eq. 01 
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Figure 4: Lowest order Feynman diagrams of lepton Savour violating 
processes induced by /121/221 couplings (see Eq. (|2. 

4.3 /i conversion in nuclei induced at tree-level 

Here, we consider the possibility that fi~ — > e~ conversion in nuclei is induced 
at tree- level. This can arise through some of the LQD terms (A'j^-^). As an 
example, we consider a model in which only X'^^i and A221 are non-zero, so 
/i~ — > e~ conversion is generated at tree-level while /i+ — > e"'"7 and /x"*" 
e~^e~e~^ are generated at one-loop level (Fig. H). 
The LFV vertices are 

rtd _ /l21'^221 aR _ /l21'^221 on^ 
^ ~ o — 2 ' "^2 — ~^T^ 2~' 1^4. ZUJ 



+J(-2-21og£|-iW)lJ. (4.21) 

As before, we set = — for fi~ — >■ e~ conversion, and = for yU+ — > 
e"'"e~e"'". The ratios of branching ratios are 

^'^ = 1.1, (4.22) 



Br(/i+ — i> e+e~e+) 

= 2 (1) X 10^ (4.23) 



R(/i- e- in Ti (Al)) ,,5 



Br(/i+ — > e+e~e+) 

Here we assume mj^ = = 300 GeV. Since /i~ — e~ conversion is 
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induced at the tree-level, its event rate is much larger than that of other 
processes, as expected. 

In /i+ — >■ e~^e~e'^, the off-shell photon penguin vertex (Af) dominates 
over the other contributions because of the ultraviolet log-enhancement. 
Therefore, the ratio of branching ratios Br(/i+ e~''7)/Br(yU+ e+e~e+) 
and the P-odd asymmetries Ap-^ and Ap^ (which are presented in Table |l|) are 
very similar to those we obtained in the previous subsection. The order one 
numerical differences come from the different sfermion masses used in both 
cases and the fact that there are quarks and not leptons running around the 
loops. Results for other similar examples are also listed in Table [1|. 

In the case of /i"*" — e^e~e~^, the fact that we choose a fixed value of q'^{= 
0) instead of integrating over all possible values leads to some uncertainty. 
These, however, are not important as far as our intentions here are concerned. 
We note that the numbers presented in Table |I] for ratios of branching ratios 
when there are first generation quarks running around the loops are uncertain 
by some tens of percent. 



4.4 Large T-odd Asymmetry in /i+ e+e e+ 

It is important to understand if any interesting effect can be obtained if 
more than a single pair of RPV couplings is present. Here we consider the 
possibility that /i+ — >• e^e~e~^ is generated at the tree-level, but that the loop- 
level contributions of on-shell (and off-shell) photons is comparable. This can 
be accomplished by having, for example, nonzero A131A231 and A133A233 3> 

Al3lA231- 

In this case, all of B^, A^, and A2 can be comparable, and there is the 
possibility that the T-odd asymmetry in yU+ e"*"e~e+ decay (Eg. ( p.8|) ) is 
large. We proceed to discuss this in more detail. 

We will consider the most general case in which all effective couplings , 
A^, and A2 are independent (as may be effectively the case if many RPV 
couplings are relevant). In this case, the T-odd asymmetry (Eg. ( p.8|) ) can be 
written as 

3(aiiCii — 0126*12) 

At ^ 



2(026*2 + 046*4 + a^C^ + ayCj + OgCg) ' 

3aiix sin(6'2 — 6^1) — 3ai2{?/ sin 62 + x sin(52 — ^1)} 



X . (4.24) 



17 



where 



X = 4a2a;^ + Aa^i^x^ + + 2xy cos 6i) + 

— 2a7Xcos(6'2 — ^^i) — 2ag{yco?,02 + xcos(^^2 — ^i)}- 

Here (02, 04, 05, ay, ag, an, ai2) = (1.8, 0.96, 88, 14, 7.5, 2.0, 1, 6), x = l^f/Af |, 
y = /4:7iaA2\, and 9i {62) is the relative phase between and (A^). 
Even when x, 7/, ^1 and 62 are treated as independent parameters, this T-odd 
asymmetry has a maximum value. 



A 



T I max 



24%, 



(4.25) 



when 



X 



01 

02 



A^ 



2.56, 



A-waA^ 
-2.28, 
-1.56. 



4.23, 



(4.26) 



This upper limit is quite general, and applies to any extension of the SM. It 
can be obtained directly from the most general effective Lagrangian which 
parametrises /x"*" — > e+e~e+ [^. 

Fig. ^ depicts the value of the T-odd asymmetry and the ratio of branching 
ratios of z^"*" e"'"7 and /x"*" — > e^e~e^, when we fix 61 = —2.28, 62 = 
— 1.56 (same as at the maximum point). As can be seen from Fig. ^ these 
two observables are strongly correlated. In the region where the T-odd 
asymmetry is relatively large, the branching ratio for /i+ — > e"'"7 tends to be 
much bigger than the one for /i"*" — > e~^e~e~^, since an on-shell photon coupling 
A2 comparable to A^ and B^ is required in order to obtain a large T-odd 
asymmetry. In this case, the branching ratio of /i"^ e^e~e^ is dominated 
by the A2 coefficient due to the relatively large collinear infrared logarithm 
(see Eq. ( 3.19| )) and we obtain a ratio of branching ratios similar to the one 
obtained for the MSSM with heavy right-handed neutrinos (Eq. ( [4.6[ )). 

In a generic RPV model, the T-odd asymmetry is unlikely to be close to 
its maximum value (Eq. (|4.25|) ) because the branching ratio for yU+ —>■ e+e~e^ 
is expected to be comparable to (an in some cases even much larger than) 
the one for /i"*" — > e+7, as we argued in the previous subsections. It is. 
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Br(|i ^ ey)/Br(|i^eee) 




8 10 



Figure 5: Constant contours of the T-odd asymmetry (top) and the ratio 
of branching ratios Br(yU+ e+7)/Br(/i+ e+e^e+) (bottom) in the {x x 

y) plane, x = \^\ and y = |^~^|- The relative phases between 

and A^, are fixed at {9i, 62) = (—2.28, —1.56). See text for details. 
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Figure 6: Constant contours of the T-odd asymmetry and the in the {m^^ x 
IA133A233/A131A231I) plane, assuming that all other RPV couplings vanishes 
and that the relative phase between the two pairs of copings is 7r/2. 



however, possible to tune the various parameters in order to achieve large 
effects. In other SUSY extensions of the SM, large T-odd asymmetries can 
also be obtained in particular regions of parameter space. For example, the 
authors of |^ discuss LFV in the case of SUSY grand unified theories, and 
find T-odd asymmetries larger than 15% in some SU{5) models. 

As an example, consider a situation where A131A213 = 10^® and A133A233 = 
while all other RPV couplings are zero, leading to At = 17% for 
: = 100 GeV. Here Br e+7) = 5 x 10"^^ 

and Br e^e~e^) = 3 x 10~^^. Note that, in order to obtain large At 

values, one is required to impose a mild hierarchy in the ratios of scalar 
masses (order 10^) and a more severe, finely-tuned hierarchy in the ratio of 
couplings (order 10^), as is illustrated in Fig. 



1.6xlO'Vi 



= 500 GeV and mg^ 



5 Conclusions 



We discussed lepton flavour violation (LFV) in rare muon processes (/i+ 
e"'"7, /i"*" — *■ e"'"e~e"'", conversion in nuclei) in SUSY models with 
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trilinear R-parity violation (RPV). Such models are interesting in the sense 
that they can accommodate neutrino masses without requiring the introduc- 
tion of extra fields to the MSSM. Natural explanations for the smallness of 
the RPV couplings have been studied [|^, and are not discussed here. 

It is well known that LFV in the charged lepton sector is a very sensitive 
probe for models with RPV, and that some of the most stringent constraints 
on RPV couplings come from LFV processes. Here, instead of concentrating 
on how RPV couplings are constrained by LFV, we study the expectations 
for LFV observables in the case nature realizes SUSY with small RPV, and 
discuss a number of different observables which may play a decisive role in 
distinguishing RPV models among themselves and from other SUSY models. 

Along these lines, we considered a number of representative cases for 
different RPV models in order to understand a number of features related to 
LFV. An important observation is that, in generic RPV models, all of the 
LFV processes considered are of the same order {i.e. the ratio of branching 
ratios is of order one), or fi~^ — > e"'"7 is very suppressed with respect to either 
/i"*" —>■ e^e~e^ and/or yU~ e~ conversion in nuclei, as is summarised in 
Table |l]. This behaviour is to be compared with R-conserving SUSY models 
with heavy right-handed neutrinos, where the branching ratio of /i^ e'^'j 
is always much larger than the branching ratio for /x"*" — > e~^e~e~^ and (in 
general) the rate for /i~ e~ conversion in nuclei. 

We also argue that the P-odd and T-odd asymmetries which can be 
measured in the case of polarised e^e~e^ decays give an extra handle 

when it comes to distinguishing different models. In particular, we discussed 
whether a large T-odd asymmetry can be generated in the case of RPV 
SUSY. 

In summary, if there is indeed low-energy SUSY with small but non- 
negligible RPV couplings, it is likely that these not only contribute to Majo- 
rana neutrino masses but also will be probed by LFV in the charged lepton 
sector. If this is the case, naively higher order QED processes, such as 
yU"^ — s> e^e~e^ or /i~ e" conversion in nuclei are at least as relevant 
as the more canonical /i"'' e^7 decay. 

Independently of what the new physics beyond the SM is, it should 
be kept in mind that improving the current experimental sensitivity of all 
LFV processes is important. We hope to discuss this important issue in a 
future publication . We conclude by stressing that there are proposals for 
improving the sensitivity to e^^ down to branchings ratios of 10~^^ 

[33| and the sensitivity to —>■ conversion in nuclei down to rates of 10~^^ 
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(see Appendix B); however, in the case of the /i"*" — > e+e~e"'", there are 
no proposals for improving the current best bound, which is aheady twelve 
years old! In view of the results discussed here, we believe that experiments 
which are sensitive to smaller branching ratios for /i"*" — *■ e^e~e^ (at least as 
sensitive as the future /x"*" — >• e"'"7 experiments) are of the utmost importance. 
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A LFV effective vertices in trilinear RPV 

In this Appendix, we present explicit expressions for the LFV effective ver- 
tices B^'^, and D^''^ in the trilinear RPV models considered in the 
body of this paper. 



A.l Photon penguin vertices 

The photon penguin vertices are defined in Eq.( |3.1| ). The effective couplings 
Af^^^ {i = 1, 2) are given by 

where ^^(^^'-^'^^ are induced by R-parity violating LLE couplings through 
a lepton-sneutrino loop and a neutrino-slepton loop, respectively. 
are generated by LQD couplings through an up-type quark-down-type squark 
loop and a down-type quark-up-type squark loop, respectively. The explicit 
expressions for the on-shell photon vertices A2 follows: 

.R{u) _ ^njKsj ji2) ( K \ /A Q^ 

^2 — T7r~2 2 '^cr 2 ' 2 ' 2 ' \^-'^ ) 
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ml. g2 



,Liu) _ .L{d) 




(A.6) 



Here the functions J^'"^^ are defined by 
jW(a,6,c) 

'^1/ /'"'(^a;^ ^ x(l-x-y) 

Jo x{l — c(l — X — y)\ + a[l — x) — by[l — x — y) 



fdyj 
Jo Jo 



x + y — cy{l — X — y) + a(l — x — y) — hxy 



When b,c 1, these functions can be approximated by 

,2 I ^3 



Jy{a,b,c) = 12(1 -a)^ ' ^ 



/ , X 1 — 6a + Sa^ + 2a'^ — 6a^ log a 

= 12(13^^ ■ 

The off-shell photon vertices Ai are expressed as follows: 



(A. 12) 



" y \ ' 



i>3 \"^i>3 "^i>3 "^£>3 



^.(.) ^ _ W^,(2) / ^ ^^.,4) 
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CLj \ eLj ^Lj VI, 

1. 



J{j k \ ajli ilRj o,Rj . 



-Jf I ^ 1 ^ (A.17) 

^ ' mj rrij rrij ' ' 

d-Rj dnj dnj 



-2.f (^.^,^)|, (A.18) 

where the functions J^'"^^ are defined by 
J("(a,6,c) 

r'd—f^, — (.+2,)(i x-y) — ^ ^ 

Jo x\l — c(l — X — yj} + a(i — — oy(l — x — y) 
Jf\a,h,c) 

, (il/ / dx , ^ , ^ — . A.20 

Jo X + y — cy{l — x — y) + a{l — x — y) — bxy 

When a,b,c <^ 1, these functions are well approximated by 

Jf(a,6,c) = -l(^l + loga + 5(a/6)), (A.21) 
J?\a,b,c) = (A.22) 

where 

_ r -I -4^^ + 2(1 + 2rf)v/r^tanh-^^, for ci < i 

- I _|_4rf + 2(l + 2d)V4rf^tan-i^, for d > i. ^'^■^"'^ 

When a^h,c and 6, c ^ 1, 
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, \ -2 + 9a - ISa^ + lla^ - 6a^loga / , 

jPiaAc) = ^^^r^^, . (A.25) 

A. 2 Tree-level vertices in the fi'^ e+e e+ process 

The tree-level vertices B^^^^ in the e+e~e"'" process are defined in 

Eq.(|3.5D. Their explicit expressions are given by 



2ml ' 

jR _ -^4*11-^*12 



2m? ^ ' 



A. 3 Tree-level vertices in the ji e conversion pro- 
cess 



The tree-level vertices for /i — e conversion were defined in Eg. ( p. 201) . 
Their explicit forms are 



= (A.28) 

dm 



= -^M}^. (A.29) 



B Constraints on R-parity- violating Couplings 
from LFV Processes and Neutrino Masses 

The experimental limits on LFV processes set tight bounds on specific com- 
binations of R-parity violating couplings. The most stringent experimental 
limit on the branching ratio of /i"*" e^e~e~^ is given by the SINDRUM 
experiment at PSI pSfiH 

Br(/i+ ^ e+e-e+)|exp. < 1-0 x lO'^^^ (B.30) 

^In order to reach the current bounds, rare muon decay experiments need to stop the 
niuons before they decay. For this reason, they are constrained to analyse fi~^ decays, since 
the /i^ is readily captured by the material present in order to stop the muons and there 
are virtually no free /i~ decays. For the same reason, one can only measure the /i — > e 
conversion rate in nuclei for the /i^. 
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The present experimental limit on the branching ratio of ^'^ e^7 process 



is set by the MEGA collaboration at LANL ||36| : 

Br(/i+->e+7)|exp.< 1.2x10-". 



[B.31) 



This limit will be significantly improved (or, perhaps, LFV will be found!) 
in the near future by a new experiment at PSI |^, which claims to be able 
to observe /x^ — > e+7 events if Br(/i"'" — > e'^'-f) > 10~^^. The present experi- 
mental bound on the conversion rate of ^~ — > e~ in ^fTi was determined by 



the SINDRUM 2 collaboration at PSI [0: 



R(/x- 



m 



Ti)| 



exp. 



< 6.1 X 10 
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[B.32) 



The future proposed (almost approved) experiment MECO claims that it 
will be able to see /i~ e~ conversion in aluminium if R(/i~ —>■ e~ in ^3A1) > 
10~^^. (More futuristic proposals claim sensitivity to values of the rate of 
/i^ e~ conversion in nuclei as low as 10~^^ ||38|| !) 

Tables |, || contain current and near future bounds on the absolute values 
of some pairs of RPV couplings, assuming that all other pairs of couplings 
vanish. 

In models with trilinear RPV, neutrino masses are generated at one- 
loop via squark (slepton) exchange for LQD (LLE) operators. Under the 
assumption that the left-right sfermion soft mass-squared mixing terms are 
diagonal in the physical basis and proportional to the associated fermion 
mass 
to IT 



mir, „ oc mfTTij:), the formula for the neutrino masses can be simplified 



fix) 



^c^ijk^ikj 

167r2 



rrii 



mi 



[B.33) 



(xlna; — x + l)/(a; — 1)^ 



Here, mf^ is the fermion mass of the zth generation inside the loop, is 

the average of the fu and /r, squark masses, and ric is a colour factor (3 for 
LQD operators and 1 for LLE operators). This expression implies that the 
heavier the fermions in the loop, the stricter the bounds For example, 
demanding rrie^ < 1 eV for sparticle masses of 300 GeV, rrih 
and TUs = 170 MeV 
drops to 2.3-10-^ 
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4.4 GeV 

< 4 ■ 10~^. For Ai22'^222 the bound 
while for "Super-Kamiokande-friendly" solutions with 



leads to A]^33A233 
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Table 2: Current (future) constraints on the R-parity violating couplings 
LLE (see Eq. ( |2. 1\ ) ) from LEV processes, assuming that only the listed pair of 
coupling is nonzero. The current (future) upper limits on the branching ratios 
are: Br(/i+ e+7) < 1.2 x 10"^^ (10"^^), Br(/i+ e+e-e+) < 1.0 x 10"^^ 
and R{fi- in Ti) < 6.1 x lO^^^ (R{fi~ e~ in Al) < \^-^^).'Rere 

we assume all the sneutrino masses degenerate with right-handed slepton 
masses, m- = 100 GeV, and we neglect left-right mixing in the charged 
slepton mass matrix. The notation (tree) indicates that the fi^ — > e+e~e+ 
process is generated at the tree-level. 



fi eee 



fi ^ e in nuclei 
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Table 3: Current (future) constraints on the R-parity violating couplings 
LQD (see Eq. ( |2. ![ ) ) from LFV processes, assuming that only the listed pair of 
coupling is nonzero. The current (future) upper limits on the branching ratios 
are: Br(/i+ ^ e+7) < 1.2 x 10^^^ (10^^^), Br(/i+ ^ e+e~e+) < 1.0 x 10"^^^ 
and R(/i' in Ti) < 6.1 x lO^^^ (R{i2- in Al) < IQ-^^). Here 

we assume all the squark masses are degenerate, with = 300 GeV. The 



notation (tree) indicates that the fx 
at the tree-level. 



e conversion process is generated 
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/i eee 
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yU ^ e in nuclei 



r' ■ ^ '■'-■^ ^^^^^^^ 
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(2 X 10-^) 
3.9 X 10-^ (tree) [|l|] 

(7 X 10-9) 
3.9 X 10-^ (tree) [g 

(7 X 10-9) 
3.6 X 10-^ (tree) [g 

(6 X 10-9) 

4.3 X 10-5 m 

(7 X 10-^) 

5.4 X 10-5 m 
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hierarchical neutrinos the bounds on certain products of RPV couphngs can 
be stricter by some orders of magnitude. 

When comparing these bounds with the ones from LFV in Tables ^ and 
^ we see that for a large number of models the bounds from stopped muon 
processes are significantly stronger than those from neutrino masses. A 
proper study of these processes therefore, can shed additional light in the 
issue of lepton number violation. 
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